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We show that the class of functions of topologically nontrivial gauge transformations in 
QCD includes a zero- mode of the Gauss law constraint. The equivalent unconstrained system 
compatible with Feynman's integral is derived in terms of topological invariant variables, 
where the zero-mode is identified with the winding number collective variable and leads to 
■ the dominance of the Wu- Yang monopolc. Physical consequences of Feynman's path integral 

in terms of the topological invariant variables are studied. 
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s I 1. Introduction 

\ s q . The anomalous decay processes in both electrodynamics and chromodynamics are described 

by the spatial integral of the product of the magnetic and electric field strength tensors. In 
non-Abelian field theory, this integral is known as the time derivative of the winding number 
functional. 

One of the main differences between QCD and QED is the fact of the non-invariance of the 
QCD winding number functional with respect to gauge transformations with stationary matrices 
(considered as maps of the coordinate space into the color group one [Q] ) . This is a result of the 
nontrivial topological structure of the group of stationary gauge transformations. The condition 
that the functional for the degree of the map is a finite number (normalization) determines the 
^ 'winding number' class of functions on which the gauge transformations act. These functions 

decrease at spatial infinity as 0(l/r). This means that the dynamic gluon fields behave similar 
to the fields of point charges (monopoles). 

This difference between the classes of functions for dynamic fields in QED and QCD is the 
principal peculiarity of the gauge theory of strong interactions. The present paper is devoted to 
a discussion of the dynamical consequences of the 'winding number' class of functions in QCD. 

The present approach is based on the fact that the winding number class of functions contains 
the zero-mode of the Gauss law constraint @, [|. This zero- mode is identified with the winding 
number collective variable M, Hi. 

We derive an equivalent unconstrained system compatible with the simplest canonical quan- 
tization scheme in the form of Feynman's path integral. From this path integral representation 
of the generating functional a unitary perturbation theory for non-Abelian gauge theories can 
be obtained (see 0, [7|, ||) which is similar intuitive Faddeev- Popov (FP) scheme ||. 

When comparing the scheme which will be developed in this paper with the instanton ap- 



proach [1C] the following basic differences are observed: (i) the approach is formulated in the 
Minkowski space instead of the Euclidean one, and (ii) the perturbation theory is formulated 
with physical states instead of the classical vacuum ones. 
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The paper is organized as follows. In Section 2 we discuss the statement of the problem in 
detail. Section 3 is devoted to the derivation of the equivalent unconstrained system with the 
zero-mode in the Yang-Mills theory. In Section 4, the generating functional for Green functions 
in the form of the Feynman integral is constructed for QCD as the basis for a solution of the 
problems of the hadronization and confinement. 

2. The winding number class of functions 

We consider the winding number functional 

2 ~ - ~ " 

AidjAk - -AiAjAk (1) 
v 

for the gluon gauge fields A^ = g^A® , in the non-Abelian SU C (3) theory with the action 
functional 

W = J d*x |i(Gg, 2 - Bf) + # 7 "(^ + 4) - , (2) 

where ip and ip are the fermionic quark fields. We use the conventional notations for the non- 
Abelian electric and magnetic fields 

Gt t = d At - D?(A)A b , Bt = e ijk (djAt + 9 -f abc A b A^ , (3) 

as well as the covariant derivative Df b (A) := 5 ab di + gf acb A1. 

The action (|2|) is invariant with respect to gauge transformations u(t, x) 

A> :=u(t,x)(A i + d i )u- l (t,x), :=u(t,x)ip . (4) 

It is well-known ||] that the fixation of the gauge in the classical equations of motion leaves the 
ambiguity in the choice of initial data for the gauge fields due to the invariance of the action 
with respect to the gauge transformations (Q) with stationary matrices u{t,x) = v(x). 

The group of the stationary gauge transformations v(x) in the coordinate space is topo- 
logically nontrivial and represents the group of three-dimensional paths lying on the three- 
dimensional space of the S'C/ c (3)-manifold with the homotopy group ir^(SU c (3)) = Z. The 
whole group of the stationary gauge transformations is split into topological classes marked by 
the integer number n (the degree of the map) which counts how many times a three-dimensional 
path turns around the 5'L r (3)-manifold when the coordinate Xi runs over the space where it is 
defined. The stationary transformations v n (x) with n = are called the small ones; and those 
with n 

A\ n) := v^X^Aiix^Xx)' 1 + L™ , L™ = v^X^diV^Xxf 1 (5) 

the large ones. 

In QCD, the winding number functional (||) is not invariant with respect to large gauge 
transformations (||) Q 

X[A^] = X[A}+M 1 [A,n]+M 2 [n} , (6) 

where 

MIA n] = ^Jd 3 x e^ k TrMAjLV] , (7) 
MM = "24^ / d ' X ^ Tr W L ? 1 ® = n - ^ 



d 3 xe 



tjk 



Tr 
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It determines the degree of a map A/2[n] (see Eqs. (3.33), (3.36) in [|TlJ). The degree of a map 
M2 [ n ] = n as the condition of the normalization means that the large transformations are given 
in the class of functions with the spatial asymptotics 0(l/r). Such a function Lf (S) is given by 

= exp(n$(x)), 4> = -iir^-f (r) , (9) 

r 

where the antisymmetric SU(3) matrices are denoted as 

X A := A 2 , \\ := A 5 , X\ := A 7 , 
and r = \x\. The function fo(r) satisfies the boundary conditions 

/o(0) = 0, / (oo) = 1 , (10) 

so that the functions disappear at spatial infinity ~ 0(l/r) but can have nonvanishing surface 
integrals in Eq. (|6|). We call the class of functions (0) the 'winding number' class of functions. 

The present paper is based on the evident fact that the dynamical field Ai and its transfor- 
mations both belong to the winding number class of functions. 

The second fact is that the winding number class of functions includes the topological ex- 
citations of the gluon system as a whole in the form of the zero-mode (Z) of the non-Abelian 
Gauss law constraint (2], || 



SA 



Df(A)G b 0l = -jS, (11) 



where = gij:^^^ is the quark current. The Gauss law takes the form of an inhomogeneous 
equation for the time-like component Aq of the gauge field 

(D\A)) ac A c = Df c (A)d A? + j% . (12) 

A general solution of this inhomogeneous equation is a sum of the solution Z a of the homogeneous 
equation 

{D 2 {A)) ab Z b = , (13) 

i.e., a zero mode of the Gauss law constraint, and a particular solution Aq of the inhomogeneous 
one (0) 

Al = Z a + ig . (14) 

It is the central point of our paper, that the zero-mode Z a at the spatial infinity can be 
represented in the form of the product of a new topological variable N(t) and a phase ^o(^) 

Z(t,x) I asymptotics = N(t)$ (x) , (15) 
where the phase &o(%) belongs to the winding number class of functions (|9|) 

$0 = -*3^/o(r) , /o(0) = 0, /o(oo) = l. (16) 

r 

It is determined in the asymptotical region by the equation 

(D 2 ) ab (^ b ^) = , (17) 
where &i(x) the asymptotics of the dynamical gluon fields A{ given in the same class of functions 

Ai(t,x) I asymptotics = ®i(x) = ~i^-^iak — fl (r) , (18) 
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with the boundary conditions 



/l(0) = 0, /i(oo) = l. (19) 

In this case, the single one-parametric variable N(t) reproduces the topological degeneracy of all 
field variables, provided the separation of the zero-mode phase factors of the topological degen- 
eracy from the topological invariant variables A 1 ,^ 1 of perturbation theory (i.e., the variables 
without topological degeneracy). This separation is fulfilled by the gauge transformations 

= U z (Z + d Q )Uz 1 , (20) 

A i = U z (A I + d i )Uz 1 , ip = U z ^, 
where the spatial asymptotics of Uz is 

t 

U z = Texp[ J ^ / i(t , ,x)]| asymptot ics = exp[iV(i)4>o(x)] . (21) 



A known example of fields (<3?j, $o) which satisfy Eq. ( |17| ) is the Bogomol'ny-Prasad-Sommerfield 
(BPS) monopole ($f PS ) with the Higgs field (® BPS ) |TJ 



*? P8 = -i&eiak*fr S <r), A BPS (r) = 1- ■ I, M , (22) 

2 r z esmh(r/e) 

r tanh(rye) r 

in the limit e — > 0, where e is the size of the monopole. For e — > the BPS monopole goes over 
into the Wu-Yang monopole [13| without singularity at the origin (as the solution of classical 



equations). 

In this case, the winding number functional (|IJ) takes the form of the sum of the invariant 
term and the new topological variable N(t) 

X[A] =X[A I ] +Mi[$i,N(t)] +M 2 [N(t)]=X[A I ]+N(t) . (24) 



This can be seen |T]| from the fact that for noninteger n = N(t) the responses of the winding 
number (||) are equal to 

M = sin[2^(t)]/2vr, 

and 

Mi[N(t)] = {N(t) - sin[27rJV(t)]/27r} , 

so that the sinus terms cancel each other in the sum (p4|). The function N(t) represents a one- 
dimensional parametrization of the transition between different maps of the homotopy group. 

The topological degeneration of initial data means that the points N(t) and N(t) — n are 
physically equivalent. Thus, the configuration space of the physical variables of QCD has the 
topology of a cylinder where the role of the degree of freedom which rotates around the cylinder 
is played by the topological variable N(t) as the zero mode of the first class constraint. This zero 
mode can be revealed by only the explicit resolving the Gauss law constraint ||. In the contrast 



to the zero mode of the second class constraint (i.e., gauge) treated as the Gribov ambiguity [14|, 
the zero mode of the first class constraint is the inexorable consequence of internal dynamics, 
which can be implicitly lost by the standard gauge-fixing scheme. 

The present paper is devoted to the derivation of an equivalent unconstrained system with 
the topological variable N(t) and to the construction of the generating functional for a unitary 
perturbation theory. 
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3. Yang-Mills theory 
3.1. Constraining 

Recall that for a long time the problem of quantization of non-Abelian constrained systems was 
considered as the greatest challenge of relativistic quantum field theory ]i"5|| . The success in the 
description of the quantum dynamics of such systems was the unitary perturbation theory in the 
form of the Faddeev- Popov functional integral used for the proof of the renormalizability of the 
Standard Model (Nobel Prize of 1999 for 't Hooft and Veltman). However, the foundation of the 
intuitive Faddeev-Popov integral for the Yang-Mills (YM) theory was achieved by Faddeev j(| 
on the way of the construction of an 'equivalent unconstrained system' compatible with the 
simplest canonical scheme of quantization in the form of the Feynman integral over independent 
physical variables 



Zy M [J] = /II { ft [d2A l^ — } ex P { iW YM + i J d 4 xJ?Af} . (25) 



[d 2 A T a d 2 El 

t,x ka=l 

The action of this equivalent unconstrained system 

W YM = J d'x [Ef -Af-\ {{Eff + B 2 k (A T ) + (d k a a ) 2 }} , (26) 
where a a satisfies the equation 

D{A)\ c dia c = -ge abc A°Ef , 

is derived by the constraining 

Wym = Wym (constraint) (27) 
the initial YM action in the first order formalism 

W YM = J d*x - \ [F* + B?] } (G a 0l = Af - D?(A)A b ) . (28) 

The constraining (|27|) means the substitution of the explicit solution of the Gauss law con- 
straint (see in details TQ) 

Df{A)F t b = (29) 



obtained by decomposing the electrical components of the field strength tensor Fq^ into transverse 
F^ T = Ef T and longitudinal F^ L = -diO a parts 



F-=Ef-d i a a ; d k E b k T = . (30) 
This decomposition is compatible with the perturbative gauge 

diAt := diAf = . (31) 

3.2. Constraining with the zero mode 

To formulate an equivalent unconstrained system for the YM theory in the winding number 
class of functions in the presence the zero mode Z b of the Gauss law constraint (f29|) 



A% = Z a + Ig; F \ = -Df(A)Z b + F \ ( {D\A)f h Z b = ) , (32) 
we shall follow two principles of the Faddeev foundation j|] : i) the constraint-shell action (|27| ) 

Wym (constraint) = W YM [Z] + Wy m [F] , (33) 
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and ii) the choice of the gauge 

Df($i)A\ := (34) 
compatible with the perturbation theory around the SU(2) monopole $>f 

hk = UzF^Uz 1 , Ai = U Z (A( + di)U^ , Af(t, x) = ^(x) + A^t, x) . (35) 

To present the action in the form (p3|), we can use the evident decomposition 

F 2 = (-DZ + F) 2 = (DZ) 2 - 2FDZ + (F) 2 = d(Z(DZ)) - 2d(ZF) + (F) 2 (36) 

and the Gauss Eqs. DF = and D 2 Z = which show that the zero mode part Wym of the 
constraint-shell action (|33] ) is the sum of two surface integrals 

W YM [Z] = J dtj d 3 x[^d i (Z a Df(A)Z b ) - di{F^Z a )\ = W° + W , (37) 

where the first one W° is the kinetic term and the second one W' describes the coupling of 
the zero-mode to the local excitations. These surface terms are determined by the asymptotics 
of the fields (Z a ,Af) at spatial infinity (15), (18) which we denoted, in the SU(3) case, by 
(N(t)^Q(x), &1;(x)). The fluctuations F^ belongs to the class of multipoles and since the surface 
integral over monopole-multipole couplings vanishes, the fluctuation part of the second term 
obviously drops out. The substitution of the solution with the asymptotics ( |T5| ) into the first 
surface term Eq. (37) leads to the zero-mode action 

W° = \ f dtN(tf , (38) 



where 

/Itt 
d 3 xd^ a (x)d^ a (x)] = (— ) 2 4vre , (39) 
9 



and 

„2 t> 



t= Ym^rlf'^) . (40) 

r v is a value of r at a boundary of a finite spatial volume. 

The action for the equivalent unconstrained system of the local excitations, 

W YM [F] = J d*x {e% • if - \ [e 2 + BliA 1 ) + [Df (d>)a 6 ] 2 }} , (41) 

is obtained by decomposing the electrical components of the field strength tensor F^ into trans- 
verse Fg? = Ef and longitudinal F$ L = -Df((f>)a b parts, so that 

F^ = E t -Df($)a b ; Df($)E b k = 0. (42) 

Here the function a b is determined from the Gauss equation ( |29| ) 

((D 2 ($)) ah + ge adc AfDf ($)) a b = -ge abc A^E^ . (43) 

3.3. Feynman path integral 

The Feynman path integral over the independent variables includes the integration over the 
topological variable N(t) 

Z^—[J] = f \{dN{t)Z YM [J u ] , (44) 
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where 

Z[J U \ = / II { II [d2Al f Ea] ) exp {iW YM + iW YM + iS[J u ]} . (45) 

t,x la=l ) 

As we have seen above, the functionals W, S are given in terms of the variables which contain 
the nonperturbative phase factors U = Uz (|2l| ) of the topological degeneration of initial data. 
These factors disappear in the action W, but not in the source 

S[J U ] = J d 4 xJ?A<>, Ai = U(Ai + di)U- 1 (46) 

what reflects the fact of the topological degeneration of the physical fields. 

Thus, instead of the instanton averaging over 'interpolations between different vacua', in 
order to remove the topological degeneracy, all Green functions should be averaged over values 
of the topological variable and all possible angles of orientation of the monopole unit vector (n = 
x/r) ( |l8| ) in the group space. This averaging leads to complete destructive interference of colored 
Green functions and corresponding color amplitudes. The complete destructive interference of 
color amplitudes can be interpreted as confinement in the spirit of the Feynman quark-parton 
duality 



4. Quantum chromodynamics 
4.1. Feynmann path integral 

QCD differs from the YM theory by the group SU(3) and quark fields. 

The decomposition of the spatial components of the gluon fields A 1 = + Af into the static 
monopole and the fluctuations can be used to define a perturbation theory with respect to 
the dynamical gluon fields Af. The adequate gauge for this perturbation theory is the covariant 
Coulomb gauge defined by 

I>r(*M = • (47) 

After the separation of the phase factors of the topological degeneracy, we obtain the similar 
path integral in terms of perturbative fields (E , A , ip 1 , tp 1 ) 

Z$—[J, r,, a] = /IJ dN(ty w °Z[J u ,r, u ,fj u ], (48) 
J t 

where 

Z[J u , V u ,fj u J]= f[lW] f[ [ ^ A ; d2Ea] )e^{W+iW + iS[J u ,r, u ,fi u )} , (49) 

t,X I 0=1 ) 

V U = UzV , (50) 

the zero-mode actions Wo and W' are defined by the expressions similar to (|37]) and ([38]) , in the 
first order formalism obtained by the Legandre transformation 

-Gpj = GoiFoi - 2 F 0i- 

The quark part of W' will be discussed later. The action for the equivalent unconstrained system 
of the local excitations 

W = Jd A x {l% -K l ~\ {El + SUA 1 ) + [Df (cD)a 6 ] 2 } + - m]^ 1 + jf Af } , (51) 
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is obtained by decomposing the electrical components of the field strength tensor F^ (simi- 
larly (|35|)) into transverse Fq? = Ef and longitudinal Ffa 1 = —Df b ((&)a b parts, so that 

FS! = E % -Df^)a h - Df(<Z>)E b = 0. (52) 

The function a b is determined from the Gauss equation ( ]TT| ) 

(D 2 (<S>)) ab + gf adc A d Dfm) ~a b = -j b otfi ; j b otfi = gf abc A^E^ + j b . (53) 



Equation ( |53"| ) can be formally solved by introducing a Green function G ab (x,y) defined as the 
solution to 

(D 2 ($)) ab G bc (x,y) = 5 ac 5 3 {x-y) . (54) 
Then a a can be expanded in a power series in g 

a b (t,x) = -J d 3 yG bc (x,m o t,o(t,y)- J d 3 yd 3 zG bc (x, y)gf cde A d (t, y)G e ^ (y, z)j( otfi (t, z) - ... 

(55) 

in analogy to the standard perturbation theory [j^j]. The calculation of the Green function 
G ab (x,y) is given in the Appendix. 

4.2. Dominance of the Wu-Yang monopole 

The action W° in the path integral (given by (p8|)) describes a free 'rotator' with the 
momentum spectrum Q 

P = NM N = (2irk + 9) , k = 0,1,2,..., (56) 

which follows from the constraint on the wave function for physically equivalent points N and 
N + 1 

*(JV + 1) = exp(i6)V(N). 
The action for the 'rotator' (p8|) compensates the action for the monopole background 



W[$ i ] = -\jdtjd*x[BZ{<S> i )f 



in the Minkowski space for finite values of the momentum of the rotator, so that the perturbation 
theory begins from the zero value of the Minkowski action, in the contrast to the instanton 
background. This can be interpreted as the dominance of the monopole background. 

For the self-dual BPS ansatz (22) the dominant value of the momentum of the rotator is 

(") 

2vr g l 

4.3. Hadronization 

In the lowest order of this perturbation theory, we can rewrite the instantaneous interaction 
term in (|5l| ) in the form of a current-current interaction 

Wmt = -\J d*x[Df(<S>)a b } 2 = \Jd 4 xJ d 3 yj^ (t, x)G ab (x, y)j b tot , (t, y) . (58) 

Following the QED perturbation theory of 'radiative corrections' one can formulate a similar 
perturbation theory of 'dynamical corrections' in QCD. This 'dynamical perturbation theory' is 
based on the decomposition of the action W for the dynamical variables ([47|) into the 'free' part 
Wo, the instantaneous interaction Wi n t ( p8| ) and the residual dynamical interaction W- in t 

W = W + W int + W mt . (59) 
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The lowest order of the 'dynamical perturbation theory' Wi n t = describes gluons and quarks 
in the monopole field and their instantaneous bound states. Next orders of this perturbation 
theory contain 'dynamical corrections' which describe matrix elements of transitions between 
these states. 

Bound states are obtained in a similar way as in QED (i.e., by the Schwinger-Dyson and 
Bethe-Salpeter equations with the instantaneous interaction G bc (x,y)) given in the Appendix. 
In the field of the monopole the instantaneous quark-quark potential is the sum of a Coulomb 
potential and a rising one. It is well-known that the latter one leads to spontaneous chiral 
symmetry breaking and to mesonic bound states. 

4.4. Frozen gluon approximation and U A (1)- problem 

The 'frozen gluon' approximation (FGA) follows from the Feynman path integral if we neglect 
the dynamical gluon fields: A? = (|47|). 

The FGA generating functional fl45|) takes the form 



Zfga[v,v] = I II {[##]} ex P + iWpGA + iS[ri, fj]} , (60) 

t,x 

where 

W FGA = W + Wmt (61) 

depends only on the quark fields. 

The bilocal linearization of the four quark interaction WW using the Hubbard-Stratonovich 



transformation leads to an effective bilocal meson action [17]. This meson action includes Abelian 
anomalies in the pseudoscalar isosinglet (r/o-meson) channel [|], ||, 18 1. 

In our case, these anomalies include the time derivative of topological variable N(t) due to 
interaction of the quark fields with the monopole and the zero-mode W' (see eq. (37)). 



Neglecting all mesonic channels except the r/o-meson one, we get in FGA nothing but the 
well-known gauge-invariant expression 

f 

' anomaly L'/j '/J ~ J ^'/V^^c 15 '/ K">"J 

where 77,77 are fermion sources, C v is a constant. In the BPS monopole selfdual field (0), ([l8|), 
(||)with G%i = NDf($)$ b , and 2ireDf ($)$o = B i($) we obtain the normalizable zero-mode 



WZnutoh V] = cj dtfj(t, 0)/ c75 7 ? (t, 0)-!-. / £>xGl v *G° , (62) 



d 3 xG a *G* = N(t) , (63) 



16vr 2 J ^ ^ 

as 

The physical anomaly term Wanomaiy m the r/o-meson channel takes the form |l 



Following Veneziano |19| , we can identify r/o with the field of the 7/0-meson at rest. The effective 
action including the anomaly term (62) and the zero- mode one (53) is then 



W eS = J dt[^ + m C v N + , (65) 

where ttt-o is the standard current quark mass contribution to the 7/0 meson mass. The diagonal- 
ization of this effective action leads to an additional mass of 7/0-meson 
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The finite contribution of the zero-mode to the mass of the ?7o-meson entails the disappearance 
of the mass of the collective topological variable and leads to a stable perturbation theory in 
the infinite volume limit, as in this limit the singularity-free BPS monopole converts into the 
Wu-Yang one without singularities at the origin. 

4.5. Relativistic covariance 

Recall that Dirac has also obtained the unconstrained form of QED in terms of gauge invariant 
variables for QED |2(J as functionals of the initial gauge fields by explicitly resolving the Gauss 
law constraint. The resulting unconstrained formulation of QED coincides with the one obtained 
in the Coulomb gauge pll] with the physical phenomena of electrostatics, 'dressed' electrons, 



and two transverse photon degrees of freedom. In QED, in terms of the Dirac variables [2C], the 
Poincare symmetry is realized which is mixing with the gauge symmetry [p^] . This mixing was 
interpreted in 1930 by Heisenberg and Pauli [^] (with reference to an unpublished note by von 
Neumann) as the transition from the Coulomb gauge with respect to the time axis in the rest 
frame (l^ = (1, 0, 0, 0)) to the Coulomb gauge with respect to the time axis in the moving frame 
= ZjJ + dilfj = (Ll )^. The Coulomb interaction has the covariant form 

W C = J dxdy^j r f(x)V c (z ± )jl(y)5(l ■ z) , (67) 

where jj = eip T fip T , z-~ = z^ — l^z ■ I) , z^ = (x — y)^ . This transformation law and the 
relativistic covariance of this formulation of QED has been proven by Zumino |2^] on the level 
of the algebra of generators of the Poincare group. 

In QED we can change variables to construct the generating functional of the Green functions 
in any gauges including the Lorentz invariant ones. The invariance of the corresponding Green 
functions under a change of variables (which generates the Ward-Taylor-Slavnov identities [||) is 
garanteed by the Dirac factors in source terms, which restore the Coulomb gauge Feynman rules 
in any Lorentz invariant gauge. So, the Coulomb interaction and electrostatics are consequences 
of the identification of the physical degrees of freedom which correspond to an explicit solution 
of the Gauss law, but not primarily to the choice of the gauge. For example, if one would omit 
the Dirac factors in the source terms in relativistically invariant Lorentz gauge formulations of 
QED, one would get the Wick-Cutkosky bound states formed by gauge propagators with light- 
cone singularities with a spectrum different Q from the observed one which corresponds to the 
instantaneous Coulomb interaction. 

In QCD, the moving Lorentz frame corresponds to the moving Wu-Yang monopole [p4| . We 
should only manage to realize this transformation law on the level of operators. 



5. Conclusion 

In this paper, we have shown that there is a collective topological excitation in the gluon spec- 
trum as the zero- mode of Gauss' constraint from the winding number class of functions of the 
large gauge transformations. 

This topological excitation leads to the dominance of the Wu-Yang monopole in the Feynman 
path integral compatible with the equivalent unconstrained system obtained in this work. 

In the field of the Wu-Yang monopole the instantaneous quark-quark potential is the sum of 
a Coulomb-type potential and a rising one. The latter one leads to spontaneous chiral symmetry 
breaking and to mesonic bound states. 

The r/o-meson mixes with the zero-mode so that after diagonalization of this low-energy 
action a mass shift of the 77o- meson is obtained which resolves the C^a(1) problem. 

1 Onc of the authors (V.P.) thanks W. Kummer who pointed out that in Ref. [^3| the difference between the 
Coulomb atom and the Wick-Cutkosky bound states in QED has been demonstrated. 
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Colored amplitudes contain additional phase factors which depend on the zero-mode. Av- 
eraging the colored amplitudes over the zero-mode parameters can lead to the phenomenon of 



complete destructive interference [|], [16 1, so that the color amplitudes disappear. The colorless 
ones of the type of the expectation values of the electroweak currents do not vanish since the 
zero-mode phase factors are absent. 



According to Heisenberg, Pauli [21] and Zumino |22j ], the relativistic covariance is established 
by a rotation of the timelike axis so that the non-Abelian Coulomb-gauge field moves together 
with the relativistic bound states. Recently, Faddeev and Niemi p3j constructed a similar 
relativistically covariant form of the Wu-Yang monopole. 

In summary, the present scheme for the introduction of topological global variables is a 
promising tool for the investigation of the challenging properties of QCD such as the meson 
spectrum, chiral symmetry breaking, quark and gluon confinement, and the Ua(X) anomaly. 
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A The Green function 

We can calculate the instantaneous the Green function 

(D 2 {b)) ab {x)G bc {x,y) = 5 ac 5 3 {x - y) . (68) 
In the presence of the Wu-Yang monopole we have 

(D 2 f b (x) = 6° b A - + ^ + 2(^d b - ^d a ) , 

and n a (x) = x a /r; r = \x\ . Let us decompose G ab into a complete set of orthogonal vectors in 
color space 

G ab (x,y) = [n a (x)n b (y)V (z) + £ e a a (x)e b a {y)Vi(z)]- (z = \x - y\) . 

a=l,2 

Substituting the latter into the first equation, we get 

cp 2d n 

-n>Vn + --rV n - -2V„ = o n = o, i. 

dz z z dz z z 
The general solution for the last equation is 

V n {z) = d n z l i + c n z l i , (69) 

where d n , c n are constants, and can be found as roots of the equation (l n ) 2 + l n = n , i.e. 



/.™ = -i±4±^; g = Zl±^±g. (70 ) 

It is easy to see that for n = we get the Coulomb-type potential do = — 1/47T, and for n = 1 
the 'golden section' potential with 

l\ = -i±5^ „ -1.618 ; £ = Zi+v? „ 0.618 . (71) 

The last potential (in the contrast with the Coulomb-type one) leads to spontateous chiral 
symmetry breaking and can be used as the potential for the 'hadronization' of quarks and 
gluons in QCD M, [Ll. 
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